Dislocation core energies and core fields from first principles 
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Ab initio calculations in bcc iron show that a (111) screw dislocation induces a short range 
dilatation field in addition to the Volterra elastic field. This core field is modeled in anisotropic elastic 
theory using force dipoles. The elastic modeling thus better reproduces the atom displacements 
observed in ab initio calculations. Including this core field in the computation of the elastic energy 
allows deriving a core energy which converges faster with the cell size, thus leading to a result which 
does not depend on the geometry of the dislocation array used for the simulation. 
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Plastic deformation in crystals is heavily related to the 
dislocation core properties |]|. As experimental investi- 
gation of the dislocation core is difficult, atomic simula- 
tions have become a common tool in dislocation theory. 
But dislocations induce a long-range elastic field and one 
has to take full account of it in the atomic modeling. This 
is even more crucial for ab initio calculations because of 
the small size of the unit cell that can be simulated. In 
this Letter, we illustrate this point for the screw disloca- 
tion in bcc Fe by showing that the commonly-used elastic 
description, i.e. the Volterra solution [H, has to be en- 
riched in order to get quantitative information from ab 
initio calculations. 

Two different methods based on ab initio calculations 
have been developed to model dislocations. In the first 
approach, a single dislocation is introduced in a unit cell 
which is periodic only along the dislocation line and with 
surfaces in the other directions. Surface atoms are dis- 
placed according to the dislocation long-range elastic field 
and can be either kept fixed or relaxed using lattice Green 
functions [2]. The main drawback of this method is that, 
in ab initio calculations, one cannot separate the energy 
contribution of the dislocation from the surface one. To 
calculate dislocation energy properties, one has to use the 
second approach which is based on full periodic boundary 
conditions As this is possible only if the total 

Burgers vector of the unit cell is zero, a dislocation dipole 
is simulated. Using elasticity theory, one can calculate 
the interaction between the two dislocations forming the 
dipole as well as with their periodic images and thus 
isolate dislocation intrinsic properties. 

We use this dipole approach to study the core proper- 
ties of (111) screw dislocations in bcc iron with ab initio 
calculations based on density functional theory using the 
SIESTA code as described in Ref. 0- The dislocations are 
positioned at the center of gravity of three neighboring 
atomic columns. Depending on the sign of the Burg- 
ers vector compared to the helicity of the original site, 
there are two different configurations, termed "easy" and 
"hard". The "hard" core configuration shifts locally the 
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FIG. 1: Screw dislocation periodic arrangements used for ab 
initio calculations: (a) T and (b) AT triangular arrangements; 
(c) quadrupolar arrangement. ^ = | [HI] for "easy" and 
I [111] for "hard" cores. A is the dipole cut vector. 



atoms such that they lie in the same {111} plane. From 
steric considerations, one thus expects this configuration 
to be less stable. The energy landscape experienced by 
the gliding dislocation is dictated by the energy difference 
between these two configurations which is a maximum for 
the Peierls barrier. It is therefore important to get a pre- 
cise knowledge of the corresponding core energies. 

We introduce the dipole in periodic unit cells corre- 
sponding to different dislocation arrays [tI . The triangu- 
lar arrangements of Figs. [IK and[lj3 preserve the 3-fold 
symmetry of the bcc lattice in the [111] direction. One 
can obtain two variants which are related by a 7r/3 rota- 
tion. We refer to them as the twinning (T) (Fig. [1^) and 
the anti-twinning (AT) triangular arrangement (Fig. [TJd) 
The last dislocation arrangement (Fig. [Tt) is equiv- 
alent to a rectangular array of quadrupoles. 

Simulation unit cells are built so that the two disloca- 
tions composing the dipole are in the same configuration, 
either "easy" or "hard" depending on the sign of the Burg- 
ers vector. Assuming that the elastic displacement field 
created by each dislocation corresponds to the Volterra 
one Q, the elastic energy stored in the simulation box 
is proportional to the square of the Burgers vector and 
therefore is the same for the "easy" and "hard" config- 
urations. The core energy difference between the two 
possible configurations is thus simply given by half the 
energy difference obtained from ab initio calculations for 
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FIG. 2: Core energy difference between the "easy" and "hard" 
core configurations. Solid symbols correspond to core energies 
obtained when only the Volterra field is considered and open 
symbols to core energies when both the Volterra and the core 
fields are taken into account (rc = 3 A). 



the same unit cell. This energy difference is shown as 
solid symbols in Fig. [2l The result depends on the cho- 
sen dislocation arrangement. According to the T trian- 
gular arrangement, the "hard" core configuration is more 
stable than the "easy" one, whereas the quadrupolar and 
the AT triangular arrangements lead to the opposite con- 
clusion. For a given arrangement, the convergence with 
the number N of atoms is proportional to N~^^^. The 
computational cost to directly deduce converged values 
from ab initio calculations is therefore out of reach. 

To understand how our simulation approach has to be 
enriched to lead to unambiguous dislocation core ener- 
gies, we examine the atom displacements created by the 
dislocation array in ab initio calculations. For all unit 
cells, atom displacements in the [111] direction, i.e. the 
screw component, correspond to dislocations having a 
symmetrical and compact core structure, in agreement 
with recent ab initio calculations in bcc Fe [J, 0, 0|. 
The screw dislocation dipoles also create displacements 
in the (111) plane, i.e. perpendicular to the screw axis 
(Fig. [3K) . Part of this edge component arises from elastic 
anisotropy. Nevertheless, when subtracting the displace- 
ments predicted by anisotropic elasticity for the periodic 
dislocation array from the ones given by ab initio cal- 
culations, one obtains a residual displacement (Fig. [3h) 
which looks like a combination of 2-dimension expansions 
centered at the dislocations. 

This is not included in the Volterra solution describing 
the dislocation elastic field. Nevertheless, going back to 
the seminal paper of Eshelby et al. , it appears that a 
dislocation can also lead to such a supplementary elastic 
field. Indeed, Eshelby et al. showed that a straight dis- 
location in an infinite elastic medium creates in a point 
defined by its cylindrical coordinates r and a displace- 



ment given by a Laurent series which leading terms are 

u(r, 9) = v\n{r) + uo{e) + ui{0)^ + O (^^^ . (1) 

Usually, only the two first terms of this series are con- 
sidered leading to the well-known Volterra solution ^lli] . 
This gives the long-range displacement induced by the 
discontinuity along the dislocation cut. 

Close to the dislocation core, the third term in Eq. [T] 
may be relevant too H, This corresponds to what 
is usually called the dislocation core field. Such a field 
arises from non-linearities in the crystal elastic behavior 
and from perturbations due to the atomic nature of the 
core. It can be modeled within anisotropic linear elas- 
ticity theory using line-force dipoles representative of an 
elliptical line source expansion located close to the dislo- 
cation core The core field is then characterized by 
the first moments My of this line-force distribution. We 
propose in the following an original approach that allows 
to directly deduce the moments My from quantities that 
can be "measured" in atomic simulations. 

In that purpose, we consider the elastic energy of a 
periodic unit cell containing a dislocation dipole defined 
by its Burgers vector b and its cut vector A. Each dis- 
location also creates a core field corresponding to the 
moments given by the second-rank tensor M. An homo- 
geneous strain can be superposed to the heterogeneous 
strain created by the dislocation dipole. This contributes 
to the elastic energy by an amount 
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where S is the area of the simulation unit cell perpendic- 
ular to the dislocation fines, h the corresponding height 
and C'ijki the elastic constants. The homogeneous stress 
is defined as 
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where the elastic compliances Siju are the inverse of the 
elastic constants. 

When the dislocations do not create any core field 
(M = 0), one recovers the fact that the elastic energy is 
minimal for an homogeneous strain equal to the plastic 
strain produced when the dislocation dipole is introduced 
in the simulation unit cell The core fields induce a 
second contribution which is proportional to the disloca- 
tion density, thus allowing to define a dislocation forma- 
tion volume. Our ab initio calculations lead for a screw 
dislocation in bcc iron to a dilatation perpendicular to 
the dislocation line, SV± = (e^^ -f e§2)5'/2 = 3.8 ±0.3 A^, 
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FIG. 3: Planar displacement map of a periodic unit cell containing a screw dislocation dipole obtained from ab initio calculations: 
(a) total displacement, (b) after subtraction of the Volterra elastic field, and (c) after subtraction of the Volterra and the core 
elastic fields. Vectors correspond to (111) in-plane displacements and have been magnified by a factor 50. Displacements 
smaller than 0.01 A are omitted. For clarity, displacements of the six atoms belonging to the cores of the two dislocations are 
not shown in (c). Atomic positions are drawn as circles with a color depending of their original (111) plane. 



and to a contraction along the dislocation line, Wj j ~ 
e%^Sl2 — —1.3 ± 0.2 A^, where the formation volumes 
are defined per unit of dislocation line. 

Instead of letting the unit cell relax its size and shape, 
one can also keep fixed the periodicity vectors and mini- 
mize the energy only with respect to the atomic positions. 
The simulation box is thus subject to an homogeneous 
stress from which the moments responsible for the dislo- 
cation core field can be deduced using Eqs.[3]and[4l The 
component a^z of the "measured" homogeneous stress is 
negligible compared to cth and 022, in agreement with 
the following argument. 

For a [111] screw dislocation in a cubic crystal, be- 
cause of the 3-fold symmetry, the tensor M is diagonal 
with Mil = and M33 = if the unit vector 63 cor- 
responds to the [111] direction. The core field is thus 
a pure dilatation in the (111) plane. This is true when 
the dislocation is in a stress-free state or if the stress ex- 
perienced by the dislocation also obeys this 3-fold sym- 
metry. The ab initio calculations indeed lead to such 
a tensor M for the two triangular arrangements. The 
quadrupolar arrangement induces a stress which does not 
obey this symmetry. Because of the moment polariz- 
ability we obtain different values for M\\ and M22 
in this case. Nevertheless, all dislocation arrangements 
used in ab initio calculations converge with the cell size 
to Mil = M22 = 650 ± 50 GPa.A^ for both "easy" and 
"hard" core configurations. As for the contraction ob- 
served along the dislocation line, it arises from the elastic 
compliance 6*1133 which couples the strain component £33 
with the force moments Mn and M22. 

Knowing the moments, we model the dislocation elas- 
tic displacement as the superposition of the Volterra and 
the core fields. We can thus compare the displacement 
given by ab initio calculations with the one predicted 
by elasticity theory for the dislocation periodic array 0|. 



Looking at the difference between the fields given by the 
two modeling techniques for the in-plane (111) compo- 
nent (Fig. [S]:), one sees that elasticity theory perfectly 
manages to reproduce the displacement given by ab ini- 
tio calculations, except for atoms which are too close to 
the dislocation cores. It is clear that the superposition 
of the core field to the Volterra solution greatly improves 
the description of the dislocation elastic field. 

The excess energy i?, i.e. the energy difference per unit 
of height between the unit cell with and without the dis- 
location dipole, is the sum of the two dislocation core 
energies E^°'''^ and of the elastic energy. 

E - 2^™'° + ~ h,K%h, In (r^) + M^jKI^^Mu \ (5) 

where if ° and are definite positive tensors which only 
depend on the elastic constants. contains the elastic 
interaction between the two dislocations composing the 
primary dipole, as well as the interaction with their pe- 
riodic images. The core fields modify this interaction 
energy as the dislocations now interact not only through 
their Volterra elastic fields but also through their core 
fields and the combination of these two elastic fields. 
The last term in Eq. [5] corresponds to the increase of 
the dislocation self elastic energy due to their core fields. 
The cutoff distance Tc is introduced because elastic fields 
are diverging due to elasticity inability to describe atom 
displacements in the dislocation core. The core energy 
that can be deduced from atomic simulations therefore 
depends on the value of r^. 

We use Eq. [5] to extract dislocation core energy from 
atomic simulations: E and M are deduced from ab initio 
calculations, whereas E'^ , K° , and are calculated with 
anisotropic elasticity theory. A core radius slightly larger 
than the Burgers vector (re = 3 A) leads to reasonable 
core energies and a good convergence with the size of the 
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FIG. 4: Sketch of the (111) in-plane displacement created by 
the triangular arrangement of dislocations in their "easy" core 
configuration. For the T variant (a), the displacements due 
to the Volterra and the core fields have the same sign and 
sum up in the region between two neighboring dislocations, 
whereas they partially cancel for the AT variant (b). 



simulation unit cell. The core energy difference between 
the "easy" and the "hard" core configurations of the screw 
dislocation in bcc iron converges now rapidly to a value 
which does not depend on the geometry of the dislocation 
arrangement (Fig. [2]). For all simulations, the "easy" core 
configuration is more stable than the "hard" one, with 
a core energy converging respectively to E^™y — 219 ± 
1 meV.A-i and E^°J^l = 227± 1 meV-A"!. 

We can now understand why the simple approach, 
where only the Volterra elastic field is considered, leads to 
core energies which strongly depend on the geometry of 
the dislocation array. Looking at the (111) in-plane dis- 
placement created by each component of the dislocation 
elastic field, the Volterra part oscillates as a function of 
9 between a compression and a tension type, whereas the 
core-field only leads to a compression. This is illustrated 
in Fig. m for the two variants of the triangular arrange- 
ment with the dislocations in their "easy" core configu- 
ration. It is clear on this figure, that the effects of the 
Volterra and the core fields will sum up in the regions 
between two neighboring dislocations for the T variant 
(Fig. UK) , whereas they will partially compensate for the 
AT variant (Fig. SJd) . One thus expects a stronger elastic 
interaction between dislocations for the T variant than 
for the AT one. This is the opposite for the "hard" core 
configuration, as changing the sign of the Burgers vector 
reverses the Volterra elastic field without modifying the 
core field. When neglecting the dislocation core field, one 
thus overestimates the elastic energy difference between 
the "easy" and "hard" core configurations for the T vari- 
ant and underestimates it for the AT one. On the other 
hand, the coupHng of the Volterra and the core elastic 
field leads to a negligible interaction between neighboring 
dislocations for the quadrupolar arrangement because of 
its centro-symmetry. This arrangement actually appears 
as the best-suited one to extract quantitative information 
from atomic simulations Q. 

This dilatation due to the dislocation core field is not 
specific to iron. When analyzing previous ab initio cal- 
culations 0, [3| we can conclude that screw dislocations 



exhibit a similar core field in other bcc metals like Mo 
and Ta. On the other hand, empirical potentials may 
fail to predict such a core dilatation. This is the case 
for Mendelev potentia l |l7 | \yhich is often used to study 
" ll. 



dislocations in iron [^, [li 



In addition to determining the formation volume of 
the screw dislocation in iron and modeling it within 
anisotropic hnear elasticity theory, our study shows that 
considering this core field is crucial when deriving from 
atomic simulations dislocation parameters like their core 
energy. This supplementary elastic field should also in- 
fiuence any energy differences, like the Peierls barriers, 
and stresses extracted from atomic simulations. More- 
over, because of the formation volume associated with 
this core field, a dislocation can interact with an hydro- 
static stress. Close to the core, it will modify the dislo- 
cation interaction with point defects (l9| . 
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